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Abstract 

We compare different formalisms for the calculation of lepton pair emis- 
sion in hadron-hadron collisions and discuss the consequences of the approx- 
imations inherent to each of them. Using a Lorentz-covariant and gauge- 
invariant formalism proposed by Lichard, we calculate lepton pair emission 
via bremsstrahlung in proton-proton and neutron-proton reactions at energies 
between 1 and 5 GeV to leading order in the virtual photon four-momentum. 
These new results suggest that some previous bremsstrahlung calculations 
based on varieties of the soft-photon approximation might have somewhat 
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overestimated the DLS low-mass dilepton cross sections. We find that the re- 
spective intensities of dilepton production through pp and np bremsstrahlung 
are energy-dependent and become comparable at E'kin ^ 2 GeV. We also 
calculate the polar angle anisotropy of the lepton spectrum. 
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I. INTRODUCTION 



Dileptons and photons are probably the best carriers of information from the hot and 
compressed nuclear matter produced in the early stages of heavy- ion collisions [|l| . In princi- 
ple, those electromagnetically interacting particles can leave the hadronic environment from 
which they are created without significant disturbances, offering a relatively clean probe of 
the nuclear collision dynamics. In heavy ion collisions at incident energies of 1-5 A GeV, 
where the Lawrence Berkeley Laboratory Dilepton Spectrometer (DLS) has already taken 
measurements the most important sources of e~^e~ pairs seem to be Dalitz and radia- 
tive decays — mainly from rj mesons and A's, pion-pion annihilations, and nucleon-nucleon 
bremsstrahlung Measurements of lepton pair production in single nucleon-nucleon 

reactions have also recently been performed in the GeV energy regime. In those, the mea- 
sured pd/pp dielectron yield ratios display a clear beam energy as well as invariant mass 
dependence [^. This suggests that the dominant mechanism for dilepton production may 
be changing as the beam energy per nucleon increases from 1 to 5 GeV. In particular, one 
should pay attention to the opening of inelastic nucleon-nucleon channels. The latter have 
been shown to play an important role for dielectron production in nucleon-nucleon colli- 
sions at 4.9 GeV To eventually understand quantitatively and completely the relative 
role of all these contributions and their excitation function in the complex environment of 
nucleus-nucleus collisions, it is vital to first calculate the lepton pair production cross sec- 
tions for individual processes as accurately as possible. Here we shall concentrate on the 
bremsstrahlung generation of lepton pairs in the case of nucleon-nucleon reactions. 

Several different calculations for electron-positron pair emission through nucleon-nucleon 
bremsstrahlung have been performed for reactions at, and slightly above, 1 GeV. Some of 
the more sophisticated approaches used relativistic one-boson exchange (OBE) Lagrangians, 
with the coupling to the electromagnetic field done by minimal substitution These 
approaches are thus entirely Lorentz-covariant and are also gauge-invariant in the electro- 
magnetic sector. Note that the OBE dilepton calculations can be made gauge-invariant even 
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when form factors at the strong interactions vertices are used 0. Such approaches have all 
the transformation properties that are required of a complete theory but they are not com- 
pletely satisfactory in two respects. First, they are very cumbersome. If several meson 
fields are involved, the number of Feynman graphs to be evaluated proliferates rapidly and 
the difficulty of the calculations increases accordingly. The coupling constants in the OBE 
model are fitted such that the total nucleon-nucleon cross sections are reproduced as closely 
as possible. This exercise thus has to be repeated for each colliding system. Second, even 
if two different OBE calculations with two slightly different set of ingredients (meson fields, 
form factors and coupling constants) can do a good job of generating total nucleon-nucleon 
cross sections that are in agreement with experimental measurements, generally they will 
have different predictions for the differential cross sections. As we shall see below, there is a 
way of writing the low invariant mass dilepton production cross section in nucleon-nucleon 
collisions such that it clearly depends on the differential elastic cross section of the colliding 
partners. This fact thus imposes very stringent requirements on the OBE models as far as 
their ability to predict lepton pair production yields is concerned. This point was recently 
made in the literature |T(]|,|5|. Because of the above considerations, the calculations pertain- 
ing to the bremsstrahlung generation of low invariant mass lepton pairs in nucleon-nucleon 
collisions have used a "soft photon approximation" ||lT| of some kind or another. Almost all 



of the recent calculations of dielectron production in nucleon-nucleon (and nucleus-nucleus) 
collisions that have used the soft photon approximation have used as their starting point a 
formula suggested by Riickl 

This equation links the cross section for production of dileptons via virtual photon 
bremsstrahlung to the bremsstrahlung cross section for real photons. In the above, p± 
is the three-momentum of the electron or positron, E± is the energy, = (p+ +p_)^ is 
the dilepton invariant mass squared, q is the photon momentum and u is its energy. The 
fine structure constant appears as a. 



The derivation of soft photon formulae in the context of bremsstrahlung emission of 
lepton pairs has recently been re-analyzed ||13|]. It was shown that Riickl's formula was 
not properly Lorentz-covariant and did not contain the relationship between dilepton cross 
section and virtual photon cross section that is required by gauge invariance |T^. In the 
context of the new analysis of dilepton production via virtual photon bremsstrahlung, a new 
set of formulae were derived for the leading order and next-to-leading order contributions to 
lepton pair emission |jl3|. 

It had been known for some time that the Riickl formula was not entirely correct llT5|,|7|-p| . 
Nevertheless, in comparison with OBE calculations, it was deemed quantitatively ade- 
quate 0,^. It is not the point of this paper to restate this known fact. We rather wish 
to investigate the empirical predictions of the new formalism in the practical framework of 
quasi-elastic nucleon-nucleon collisions at energies around and above 1 GeV, to establish 
whether the new leading term approximation deviates significantly from previous works. 
We also plan to investigate quantitatively and systematically the effects of different layers 
of approximation found in several calculations. Our paper is organized as follows: in section 
|T| we introduce the general formalism, in section ^ we compare different approaches to the 
generation of soft virtual photons and discuss the differences. In section |V| we explore the 
consequences of a better treatment of the many-body phase space. In section 0, we discuss 
how some of the approximations seen in sections |I| and |T| can be put into perspective by 
an angular anisotropy analysis of the lepton spectrum. We then summarize. 



II. GENERAL FORMALISM AND HADRONIC ELECTROMAGNETIC 

CURRENTS 

Consider the reaction 

a + b ^ c + d + e+e" (2) 

where a, b, c, and d represent nucleons. The generic Feynman diagram for the leading- 
term contributions to the emission of a soft virtual photon are shown in Fig. |l|. This is 
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the process on which we concentrate in this work. See the appendix for a mathematical 
description of what is a leading term and what isn't. The circle schematically represents the 
strong interaction. In the limit of soft photons, real or virtual, the radiation from the strong 
interaction blob is a sub-leading contribution. When a, b, c and d are all protons, there are 
eight Feynman graphs to be added coherently for the calculation of lepton pair emission. 
Note in passing that there have been arguments that bremsstrahlung from pn reactions 
should be significantly more important than that from pp. These were based on the fact 
that, nonrelativistically, the first non-vanishing multipole contribution for pp appears at the 
quadrupole level, whereas for np it is at the dipole stage. It has however been recently shown 
that such arguments do not hold for relativistic collisions 0. In fact, the lepton pair yields 
from these two processes are comparable at 4.9 GeV 0. There will be more on this later. 
Naturally, the pp bremsstrahlung contribution is crucial in the interpretation of the ratio of 
dileptons produced in pd and pp reactions. We will discuss pp contributions in this work, 
also. 

We evaluate and sum the amplitudes of the relevant Feynman diagrams of the type in 
Fig. |I|. Anticipating the soft photon limit, we assume that the hadronic part of the total 
matrix element is unaffected by the fact that one of its legs is slightly off-shell. We write 
this on-shell matrix element for elastic nucleon-nucleon scattering as A^o- We omit the 
momentum labelling of the initial and final state, for simplicity. Squaring the net matrix 
element, A4, and summing over final state spins and averaging over initial states, we obtain 
to leading order (see appendix) 
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In the above, 



^ J2 \Mf = ATTalMofj^ri,, . (3) 



is the hadron electromagnetic current, and the lepton tensor is 

L'" = ^(2b+V+p_>+'^)-MV^) . (5) 
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The Q's and p's represent charges and four-momenta for the particles, and = {p-^ +p- 



Note that J^q^ = 0, as a consequence of gauge invariance. |A^oP is the on-shell hadronic 
elastic scattering matrix element, squared, summed over final spins and averaged over initial 
spins. 

After performing the appropriate contractions, we may write the differential cross section 
for e~^e~ pair production with invariant mass M and energy go as 

d^p^d^p. ~ mT^M^ R,{s,ml,ml) i ^""^ ~ M^^^ ' \^^d^cdt , (6) 
where (p* is the azimuthal angle of one of the outgoing hadrons in the center-of-mass frame 
and the four- vector I = p^—p^ is the difference of positron and electron four-momenta and t 
is the four-momentum transfer. Equation follows from Eq. (2.15) of ||13| after specifying 
the nonradiative cross section as 

dao = -^dn, = ^^^d^pjt , (7) 



and including the phase-space correction introduced in [|T^. 

In the evaluation of the original Feynman diagrams, we have neglected the four- 
momentum q of the virtual photon in the phase-space 6 function in order to reproduce 
the kinematics associated with the on-shell elastic differential nucleon-nucleon cross sec- 
tion, da/dt ||Tl|. Because of this approximation, we include fl^ the ratio of two-body phase 



space R2{s2, ml, ml) / R2{s, ml, ml) evaluated at S2 and s, where s is the invariant energy 



squared available for all the final-state particles and S2 = s+M"^ — 2qo*^/s. This approximate 
correction prevents the lepton pair from violating the overall energy-momentum conserva- 
tion laws and thus has a significant effect on dilepton distributions. The ratio constructed 
from Eq. with and without the correction factor correctly drops monotonically to zero in 
the limit of maximum invariant mass. Handling the phase-space properly is quite important. 
However, setting g = in the four-dimensional delta function also forces the expression for 
the current to be evaluated on-shell, i.e. Pc + Pd = Pa + Pb- Clearly, this will have some 
effect on the value of the current of Eq. (Q). For the purpose of clarity, we will investigate 
separately the effects of a complete treatment of phase space in section 



An exercise in relativistic kinematics gets one from Eq. (j^) to the Lorentz-covariant 
differential cross section in global dilepton variables 



where /i is the electron mass. Our usage of Riickl's approach for dilepton production consists 
of neglecting the term {p+p'^ + P+P-) of the leptonic tensor p|, Eq. (|^), and analytically 
continuing the real photon energy to go = Vo^ + M'^ ■ This procedure is however not Lorentz- 



covariant [|T^]. One then obtains Eq. (§) multiplied by an additional factor of 3/2 ( 1 + 
2/i^/M^)~^. From here on, we shall discuss electron-positron pair production exclusively 
and we systematically will use /i = 0. 

To complete our analysis of previous efforts we now turn our attention to the electro- 
magnetic current. For emission of real photons in reaction the current is 

= -Q^^ - Q^JPL + + Qd^. (9) 

Pa- q Pb- q Pc- q Pd- q 

This approximation was also used by several authors, in the case of soft virtual photons 
r5|,|3|JT8|J^j5[1 . It is worth mentioning here that for the case of virtual photons, the "soft" 



limit is not unambiguously defined: for virtual photons, what is meant by a "soft photon 
limit"? Here, we use this term in connection with the condition M — 0. However, bear in 
mind that the lepton pair energy and lepton pair three-momentum can individually still be 
quite large. Thus M — > is in fact not a sufficient condition to neglect both go and q in the 
delta functions, even if this omission can partially he corrected [e.g. Eq. (P)). The opposite 
limit where M is large implies that the lepton pair energy is large and most certainly can't 
be ignored in the phase space delta function. By the same token, off-shell effects on the 
strong interaction amplitude can be small in the M — > limit, but won't be, in the large M 
region. In this work we do not include a discussion of off-shell effects. We should however 
keep this point in mind when discussing the regions of validity of our calculations. 

There are thus several possible ingredients that can affect the results of lepton pair 
emission calculations of the "soft photon type" . One resides in the use of the Riickl formula, 
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a widely-used expression which nevertheless is not Lorentz-covariant. This first point is 
easily settled: our interpretation of Riickl's formula for electron-positron pair production 
instead of Eq. (^) simply introduces an overall factor of 3/2. Another possible source of 
disagreement between calculations has to do with the use of the current of Eq. (|^) as the 
electromagnetic current for hadrons when massive lepton pairs are emitted. For lepton pairs 
with very small invariant masses the current represented by Eq. (0) of course reduces to 
that of Eq. (P). Also, the exact treatment of dynamics will play a role. For the time being, 
we set g = in the four- dimensional Dirac delta function and simply use the phase space 
correction factor associated with Eq. (^. In the next section, we investigate the effects of 
using different hadronic electromagnetic currents. 



III. COMPARING CURRENTS 

Squaring the current (^, as is required by Eq. (|), we obtain an expression that depends 
on the spatial orientation of the virtual photon. Taking an angular average, and for an equal- 
mass reaction of the type specified in Eq. with rua = nib = rric = = we arrive 
at 



-J' 



-2XlQM^ 

s s 

-2XlQM^ 

s s 

+2XlX2{QaQc + QbQd) 

S 

+2X,X2{QaQd + QbQc)^ ^F{^i , 



(10) 



where the function 



F{x,y] 



2VR 



In 



x-y-y 



R] 



[x ■ y — x"^ + ^/R] [x ■ y — y^ + a/R] 



(11) 



and where the radicand is 
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R = {l-x-yf-{l-x'){l-y^). (12) 

To make the formula more compact, we have introduced the variables 

1 1 

The velocities (3 are related to the above definitions and to invariants through 

/3:^=/?;^ = Ah(i-^), /?:^=/3:^ = Ah(i-^), 
/3:-/3;' = -A?(i-^), ;3:./3; = -A^(i-i!^), 

s s 

z^: ■ = -z^: ■ = aia^tii - ^ + -). (m) 

s s 

Notice that terms directly proportional to 1/go have exactly cancelled and disappeared upon 
squaring the current. If we set M = in Eq. (|10|), we recover Eq.(3.6) of Ref. [|T9| . 



We have calculated the ratio Rq = — ^virtual / " ^reai with the angular-averaged currents 
of Eq. (^) and Eq. (^. This ratio depends on the invariant energy s, the four-momentum 
transfer t, and the dilepton invariant mass and energy: M and go- For the purpose of com- 
parison, we arbitrarily input some reasonable values and plot Ro against the virtual photon 
three-momentum in Fig. ^. We find sizeable deviations from 1 in this quantity. The ratio 
for pn processes is sensitive to the value of ^^c.m.; the scattering angle in the nucleon-nucleon 
centre of mass, whereas it is quite insensitive for the pp case. Interference effects evidently 
play a role here. To better understand the effect of the "complete" current (also labelled 
"virtual" , above) on the differential cross section for dilepton production, we integrate over t 
and plot Ri = (^]^^3^) virtual/ (^]g^)reai as a function of the cm. virtual photon momentum 
|g I in Fig. ^, for different dilepton invariant masses. We are now essentially comparing the 
effect of using Eq. (^ over Eq. (^ in a dilepton production calculation. Also, to actually 
compare two sets of calculations that have been done previously, the numerator is derived 
from Eq. (^, while the denominator is derived from Riickl's approach. 

Before discussing the results in Fig. |^, let us restate that the net lepton pair spectrum 
will depend on the details of the differential nucleon-nucleon elastic cross-section, da/dt. 
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This is evident from Eq. (j^) and Eq. (^. It thus follows that the present analysis will have 
to be redone for different systems at different energies, with possibly different conclusions. 
For example, quark-quark and pion-pion virtual bremsstrahlung calculations have attracted 
some recent attention in connection with ultrarelativistic heavy- ion collisions |T9| , |20[] . The 
consequences of the present formalism on such studies is presently being analyzed ^1 



It is thus imperative to have an accurate parametrization for the elastic nucleon-nucleon 
differential cross section da/dt over the relevant range in t. For pn elastic collisions, da/dt 
is nearly symmetric for kinetic energy less than 1 GeV. But at higher energies, the observed 
distributions are not symmetric about Oc.m. = 90° but rather develop a stronger forward peak. 
This asymmetry increases with the scattering energy and can suppress the pn bremsstrahlung 
contribution to dilepton production by a factor of 4 at 4.9 GeV |llO|,^. A fairly detailed 



parametrization at 4.9 GeV kinetic energy was used in Ref. to calculate the absolute 
sizes of the pn and pp bremsstrahlung contributions. For practical reasons, we instead 
adopt the parametrizations of Refs. and for the pn and pp elastic differential cross 
sections. These two functional forms can fit the experimental data for energies up to 6 GeV 
with the necessary accuracy. 

Analyzing Fig. ^ the ratio Ri remains 2/3 for small invariant mass M < 10 MeV. 
For higher invariant masses, Ri is greater than 2/3 for almost all values of the virtual 
photon momentum, in the case of pn scattering. Increasing M in the pp case, Ri remains 
still considerably smaller than 2/3 at low momenta and drastically increases as the three- 
momentum grows. The low-momentum value of -Ri is thus quite different in the two reactions 
at hand. It is not simple to draw clear physical conclusions for these behaviours, but the main 
point is that these features suggest that using the "complete" current for virtual photons 
changes the distribution of the dielectrons in phase space, especially in the pp case. 

Upon integration, we arrive at differential cross sections in terms of lepton pair invariant 
mass. In Fig. ^ we show two curves: the dashed curve corresponds to using the Riickl 
formula with the real photon current, Eq. (^. The solid curve corresponds to using Eq. 
(1) with the current of Eq. (^). In each panel three kinetic energies are shown: 1.0, 

11 



3.0 and 4.9 GeV. We are thus in effect evaluating the difference between two approaches 
to calculating the dielectron production in nucleon-nucleon collisions: the one based on 
Riickl's equation and the real photon current (we remind the reader that this combination 
has become well known and widely used), and the one based on the more recent Lorentz- 
covariant formalism with a more general electromagnetic current for hadrons. At these 
energies (which correspond to DLS measurements) one realizes that the Lorentz-covariant 
and gauge-invariant formulae that we adopted in this work (solid curves) reduce the dilepton 
production cross section only slightly for M < 0.5 GeV in both pn and pp processes compared 
with previous results. For M > 0.5 GeV, the situation is a little more complicated as a slight 
decrease or a slight increase is observed, depending on the beam energy. We have verified 
that the ratio ( ) virtual /(^^) real is equal to 2/3 in the small M region for both pn and 
pp cases, and found it to increases monotonically with M for pn process, but decreases 
monotonically for pp as one goes toward the kinematical limit. Recall that other physical 
processes dominate the dielectron yield for M > 0.5 GeV at large enough beam energy ||^. In 
the large invariant mass region. Fig. 4 suggests that the previous bremsstrahlung calculations 
underestimate the dilepton yields at most by a factor of 2 even at the kinetic energy 4.9 
GeV, even though our formalism is not strictly applicable there. For pp scattering, it is 
worthwhile to point out that the radiation intensity is strongly depleted due to interference 
effects near the maximum M. 

Thus our findings at this point are essentially these: the "improved" formalism does not 
affect very much the soft lepton yields previously calculated. We do not address in detail the 
issue of hard leptons (M > 0.5 GeV) because there, the present formalism is inadequate and 
also because the signal is dominated by other sources. In the small invariant mass region 
M < 0.5 GeV, we find that this new formalism only amounts to the numerical difference of 
2/3. 

We have made above a formal comparison of two approaches. However, we did discover 
that the two formalisms would produce different phase space distributions of lepton pairs: 
see Fig. Let us now consider the specific case of the DLS data. It is a known fact that 
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any theoretical calculation attempting to reproduce the DLS data should first be filtered 
by the DLS experimental acceptance. It is then conceivable that the predictions of the two 
approaches could be affected quite differently by the experimental filter. We have run our 
two sets of calculations through the DLS experimental filter |^3[. In Fig. |^ we show the 
invariant mass dependence of the absolute differential cross sections, at kinetic energies 1.0, 
3.0 and 4.9 GeV. This figure now suggests that the previous calculations overestimated the 
low-mass dielectron yields by a factor of 2-3 for np and a factor of 2-5 for pp bremsstrahlung 
in collisions at 4.9 GeV. Thus the effect of the DLS acceptance is to make the differences 
between the two approaches more pronounced. The effect is smaller for lower kinetic energies. 
Definite conclusions for comparison with the data are presently being drawn as many-body 
bremsstrahlung has been shown to be quantitatively important at 4.9 GeV 0. Even though 
we feel that a complete reassessment of previous results is probably not warranted, surprises 
can not be excluded. 

Another useful comparison is the relative radiative intensities of pp and np scattering. 
As mentioned above, pp bremsstrahlung has often been neglected because of a classical 
multipole argument. We calculate the ratio here, with the formalism based on the leading- 
term approximation, i. e. with Eq. (j^) and Eq. (^. The relative intensities depend on kinetic 
energy as shown in Fig. ^ which presents the ratio R = (^^^)pp/(^^)pn as a function of 
invariant mass at different kinetic energies. Proton-proton bremsstrahlung becomes more 
and more important as the kinetic energy increases. At 2 GeV, for instance, R is nearly 0.5 
for small masses whereas at 4.9 GeV it becomes larger than 1. We shall return to this ratio 
in the next section. 

IV. ONE STEP FURTHER: THE EXACT TREATMENT OF PHASE SPACE 

In all of the above comparisons, we have consistently set g = in the 4-dimensional 
phase space delta function. We now avoid making this approximation and investigate the 
consequences. The many-body Lorentz-invariant phase space can be expressed in terms of 
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Mandelstam-type invariants [|T^, or in this case one can perform the integrals directly p4[ . 



Going back one step, let's write the general equation that leads to Eq. (^), without its phase 
space correction. In the process a + 6 — > c + d + e^e~ , the differential cross section for lepton 
pair production with invariant mass M and energy go is 

_ 1 1 r , 1 , 

xW(2.)V(p„+p.-p.-p,-,)p^p|^. (16) 



where |A1o|^ is the on-shell matrix element for the scattering a + h — > c+d, squared, summed 
over final spins and averaged over initial spins. The differential cross section can be rewritten 

as 



^° dM^d?q 4:EaEi,\wa - Vfe| 487r5 M^^ M^l 

X /(-J^)W5^(p.+p.-Pe-p.-g)^^, (16) 

where /i is the rest mass of an individual lepton. In the a + b centre of mass frame, one may 
perform some of the integrals and get 

MPa + P. - - P. - g) ^ ^ = ^ y 5 [cos e,,, j 

X ^dEddcosep^qd(j)p^ . (17) 

Furthermore, in the on-shell limit, the squared matrix element summed over final spins and 
averaged over initial spins can be related to the on-shell differential elastic cross section by 

= 167rs(s - 4m') ^ • (18) 

Integrating further, one obtains 



x^(co=.„,-i±^ll±l^f|j^^!<^^) 

dqodEddQgd COS 9p^gd(f)p^ . (19) 
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The range of the integration variables is such that the condition 

s + M^ + 2Edqo - 2y^{Ed + go) 



cos 9p J 



< 1 (20) 



2bd||g1 

has to be satisfied. 

With the help of gauge invariance, it is not very difficult to write the 4-current squared 
in terms of ten scalar products pi ■ pj. Setting = nib = = = rn, we obtain 



{2pa-q-M^y {2pb-q-M^y {2p, ■ q + M^Y {2pd ■ q + M^y 

2QaQb{^Pa ■ Pb - ^ 2QaQc{^Pa " Pc + M^) 



+ 



i2pa ■ q - M^){2pb -q-M^) {2pa ■ q - M^){2p, ■ q + M^) 

2QaQa{Apa-pd + M^) 2QbQMpb-Pc + M^) 

{2pa ■ q - M^){2pd -q + M^) {2pb ■ q - M^){2p, ■ q + M^) 

2QbQd{Apb ■ Pd + M^) 2Q,Qd{Ap, ■ p^ - M^) 

{2pb ■ q - M^){2pd -q + M^) (2p, ■ q + M^)i2pd ■ q + M^) ' ^ ^ 

In the a + b centre of mass, the scalar products can be written in terms of the integration 



variables of Eq. (p!9|). For instance, writing \pa\ = \pb\ = \p\ = ys/A — m? and the energy 
Ea = Eb = E = v/s/2, then pa ■ Pb = E^ + jpj^; Pa ■ q = Eqo - \p\\q\ cosOp^g-, Pd ■ q = 
Edqo - \pd\q\ cos 9p^g; and Pa-Pd = EEd - \p\ \pd\ (cos 9p^q cos 9p^g + sin 9p^g sin 9p^g sin (ppj. All 
other scalar products involve the above four. 

On Fig. 1^ we display three different curves for the cases of pn and pp scattering, at 
energies of 1.0, 3.0 and 4.9 GeV. We compare calculations done with Eq. (|19D, calculations 
done using Riickl's formula with the current of Eq. (^, and calculations done using Eq. (H) 
together with the current of Eq. (|). The three approaches display behaviors that are quite 
similar and are not really remarkably different. The largest deviations occur at the lowest 
bombarding energies at low M. Recall that the "soft photon limit" is not really properly 
defined by M — > 0. Running our results through the DLS acceptance filter, in Fig. ^, one 
realizes that the results obtained with Eq. ([I9|), and those obtained with Eq. (||) together 
with the current of Eq. (^ are very similar for the pn reactions. The only exception there 
is the low energy pn spectrum at low invariant masses. In the case of pp, the filter effects 
are such that the results with the complete phase space are lower than the previous ones at 
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high invariant masses. The effect of the DLS acceptance brings us to the conclusion that 
the position of the lepton pairs in phase space is mostly determined by the nature of the 
electromagnetic current, once the on-shell nucleon-nucleon cross section is fixed. The pp/np 
ratio does however show some sensitivity to the treatment of phase space, as displayed in 
Fig. 0. 

There is no question that the current of Eq. (^) is more general and appropriate for 
leptons than that of Eq. @. In a given reaction the differences in dilepton invariant cross 
sections, before filtering, that follow the use of those two currents are not minuscule but can 
not be called spectacular. There is an observable that has a remarkably different behavior, 
whether or not the current is used in calculations. We elaborate more on this in the 
following section. 



V. THE ANGULAR DISTRIBUTION OF BREMSSTRAHLUNG DILEPTONS 

We have already mentioned that the nature of the electromagnetic current for the 
bremsstrahlung producing hadrons will greatly affect the position of the leptons in phase 
space. This already became clear in our analysis of Fig. ^ Thus, it is very reasonable 
that the differences in currents could be highlighted in a treatment focusing on the fine 
points of angular distributions, for example. Angular anisotropics have recently been put 
forward as a means of distinguishing between competing lepton pair production sources ||25|| . 
This argument has power only if the angular distributions of those sources can reliably be 
calculated. Let's elaborate on this below. 

Owing to collision dynamics, the polarization of the virtual photon eventually converting 
into a lepton pair may be such that, in the rest frame of the dilepton, the single lepton 
distribution may not be isotropic. This is the essence of the idea. Following p5|, we write 



the differential cross section for emission of a lepton pair of invariant mass M, with a lepton 
coming out at a polar angle 9 in the rest frame of the lepton pair as 
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This enables us to write the polar anisotropy coefficient, B, as 

_ S{M,9 = 0^) _ 

^ ~ s{M, e = 90°) ^ • ^^^^ 

Since the full phase is obviously crucial to the proper kinematics, we use the methods of 
the previous section. The cross section to be used in the definition of 5*, above, is obtained 
by noting that for electron-positron pair production [|l3i . 



E E ^^^^^^ (24) 

where go is the lepton pair energy and i7+ is the solid angle element for positron emission in 
the rest frame of the lepton pair. With the help of the above equation, one can numerically 
integrate Eq. ([T5|) into the required format. In the on-shell approximation we are using for 
the strong matrix element, where M.q has no dependence on g, the anisotropy coefficient 
B only depends on kinematics and on the electromagnetic current. Details of the strong 
interaction do not influence it. We have numerically conffimed this by varying the strong 
interaction differential cross section and observing the constancy of B. 

We ffist investigate the polar anisotropy of the lepton distribution using the current 
associated with Eq. (|^). The results are shown in Fig. |TU|. We plot the coefficient S as a 
function of invariant mass, for different incident kinetic energies. Those energies are 1, 2, 3, 
and 5 GeV. They can be readily identified in the figure by their kinematical limit moving to 
higher invariant mass as the beam energy grows. The calculations are for bremsstrahlung 
in pn collisions. Our results at 1 and 2 GeV reproduce those of Ref. [^. Repeating the 
calculation with the full current of Eq. (H), we see our results change drastically, both 



qualitatively and quantitatively. This is now plotted in Fig. |TT[ The anisotropy coefficients 
no longer cross zero roughly together at about M = 0.25 GeV. The crossing points have 
spread out, and an important feature is that the curves have all been shifted upwards. The 
anisotropy coefficients are now mostly positive, except in the high invariant mass region of 
the last two incident kinetic energies. Also, the minima have been shifted to higher invariant 
masses. Thus, it is clear that any physical interpretation relying on the angular anisotropy 
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of the lepton spectrum will thus depend on the details of the calculation. In this observable, 
the differences in electromagnetic currents clearly stand out. 

For the case of pp bremsstrahlung, the use of the complete current is absolutely crucial. 



The current of Eq. (||) is used to compute the polar anisotropy coefficient plotted in Fig. [12. 
The current of Eq. (Q) is used to compute the curves appearing in Fig. The difference 
is striking. With the complete current, the pp signal bears some resemblance to its pn 
counterpart, whereas the current associated with real photons yields a complete different 
picture, both qualitatively and quantitatively. 



VI. SUMMARY 

In summary, we have used a Lorentz-covariant and gauge-invariant formalism to calculate 
lepton pair production cross sections via pn and pp bremsstrahlung, in the leading term 
approximation. The exact squared current derived for virtual photons is suitable for any 
energy and for any two-body equal-mass system. With the appropriate generalizations one 
could apply it to an unequal mass system and to reactions involving many-body final states. 
This is important and should be done. We find that using the complete current for virtual 
photons will produce differential dilepton cross sections that are not very different from 
those obtained using the real photon current. Since the exact virtual photon current changes 
the momentum distributions of the dileptons, the differences are somewhat accentuated by 
running the calculations through the DLS acceptance filter. However, the different phase 
space population associated with the different currents emerge in a much more striking 
fashion when plotting the polar anisotropy coefficient for leptons in the lepton pair rest 
frame. We feel that the lepton angular anisotropy might in fact be useful in the identification 
of different sources, but only once the calculations are deemed to possess the required level 
of reliability and precision. 

We do not suggest that the DLS bremsstrahlung analyses have to be completely re- 
done. However, larger differences in the approaches (the one based on Riickl's formula and 
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the Lorentz-covariant one) do show up in the pp channel, comparing with the pn channel. 
Thus to pin down the various mechanisms for dilepton production more precisely and to 
compare quantitatively to data, it does seem necessary to scrutinize some of the previous 
bremsstrahlung estimates, as our results suggest that some of them might have somewhat 
overestimated the low-mass yields. 

The relative bremsstrahlung intensities from pp and pn scattering have been an actively 
debated issue for some time. We provided a practical and quantitative way to illustrate 
this, for kinetic energies up to 6 GeV. Our results suggest that pp bremsstrahlung is quite 
important when the kinetic energy is higher than roughly 2 GeV. 
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APPENDIX: DILEPTON RADIATION FROM FOUR-FERMION 

INTERACTIONS 

In this appendix we derive the matrix element for lepton pair radiation for a generic form 
of the strong four-fermion interaction, where we take all the fermions to be charged for the 



sake of generality. We shall use and generalize slightly the methods of Ref. |]T3[. Our results 
are in agreement with those of that reference. 

We thus start with 4 diagrams like that of Fig. |1| where we consider radiation from the 
external lines, in both the initial and final state. Extra diagrams needed to antisymmetrize 
the final state can be added without problems. A little more specifically, we consider a 
process 
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a + b^c + d + r + r (Al) 

where a,h,c,d&xe all charged spin-1/2 fermions, and Pa+Pb — Pc+Pd + Q with q — p+ +p-. 
r represents the generic hadronic interaction together with external line factors {u^^\p) or 
u^^\p)) for one incoming and one outgoing fermion. For example, Tea includes the factors 
■u(*=)(Pc) and u'^^''^{pa), whereas T^b includes the factors u^^'^\pd) and u'^^'''> (pb) . If we consider 
bremsstrahlung radiation from particle a, the Feynman amplitude for that process can be 
written as 

Ma = u^'-\p-){-z{-e)j,)v^^^\p+) (^^j (A2) 

XU^'^\p,)T,,{pa-q) (:^7^^-^^p^) {-tQae7i.)u^'''\Pa) • 

Generally, F depends on a four values of incoming and outgoing four-momenta. We have 
made explicit its dependence on the four- momentum involved with the emission of the virtual 
photon, only. We can write 

Ma = -QaeL,u^"^\pc)rab{Pa - q) 2pJ- g - M^^^" ) ' (^3) 

where 

L,^^u^'-\p-h,v^'^\p+). (A4) 
Similarly, for bremsstrahlung radiation from particle b,c and d, 

Mb = -QbeL,u(''^\pd)ra.(Pb - q) { j2p^~q-M^' ) ^'""^^'^^^'^ ' ^^^^ 
Mc = QeeL,u(^^\pe)r (^^^^y^f^) r<i;>(Pc + q)u^'^\pa) , (A6) 

and 

Ma = QaeL.u^'^^ {pa)Y ' ^^^tlp' ) ^^"^^'^ + ^^""^'^^ ^^"^ ' ^^^^ 
Examine Ma- Using the Dirac equation with gamma-matrix algebra, one can show that 
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Hence, 



[7^,7 ■ q] 



2{2pa-q-M^)J 

xu^''^\pa) . (A9) 

Similar expressions are derived for the other matrix elements. The matrix element for 
reaction a + b^c + d + l^ + l~ is 

M=Ma + Mu + Mc + Md . (AlO) 

In order to identify the terms order-by-order in q, one needs to go to next to leading order. 
This is because of the fact that in the case of radiating fermions, the radiations from the 
external legs contain the leading order part and also contain some next-to-leading order 
contributions []T1| . One then systematically has to go to the next order to collect the missing 
pieces of the amplitude that will restore gauge invariance. To next-to-leading order one can 
then write 

rfe±g) =r±g"^. (All) 

Also, in the spirit of the above discussion, we introduce the electromagnetic contact term 
associated with the i^^ fermion line leaving the strong interaction core 

Ct = -eQ.^ . (A12) 

Then we have, up to leading and next-to-leading order in q, 

+u^''\p,)K:!y'^\p,) (A13) 

with 

U^'^\p,)r,,U^'^\p,) = U<^'''\p,)r,,U^'^\p,) = Mo{Sa,Sb,S,,Sd) . (Al4) 
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In Eq. ( |A13|) , the current is defined as in Eq. (|), and A^ol^a, Sfe, Sc, s^) is the matrix 
element for the hadronic on-shell scattering where Pa{.Sa) + Vbisb) Pc{sc) +Pd{sd)- The 
momentum labelhng of Ado has been suppressed for simphcity. Also, 



4pa ■ q 



-Qa {-g'^-q^ + g'^q 



Apc ■ q 

^ Pa,a dV^ 



db 



-Qc (g'^-q^-g'^q 



Pa ■ q dpi 

a\ Pc,a d^db 



Pc ■ q dp 



(A15) 



and 



ca 



Apb-q ' ' ' Apd ■ q 



+Qd {g^"q^ - g^^q' 



Pb ■ q dp^ 

Pd,a dTca 

Pd -q dpi 



(A16) 



Squaring the net matrix element, we sum over the final state spins and average over the 
initial state ones and obtain 



1 

4 



\M\\ 



where 



(A17) 



(A18) 



+ (/i < > u) 



(A19) 



= -Ti [{'yp, + m,)Tdb{l- Pa + ma) K%] , 



(A20) 



= \TT[{^-p, + m,)K^,,{^-pa + ma)T'J , 



(A21) 



and 
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B"" = |Tr [(7 ■ pd + rud) T^a {l ■Pb + rrib) 



'1 

caJ 



We have further defined 



and 



where 



ca 



Pa,a dV^f^ 



-Qa {-g'^q^ + g'^q'' 

+Qc{g'y^-g . y „ ^ , 
^ ^ Pc-q dp'i 



Pa ■ q dp^ 

a\ Pc,a dT'df, 



b-q,r]r'^ + -S^J'cab-q,r] 



^ ^ Ph-q dpi 

^ ^ Pd-q dp'l 



(A22) 



(A23) 



(A24) 



(A25) 



Next, combining the traces, one can write 

Tr [(7 ■ + Tdb {l ■Pa + ma) K'^^] + Tr [(7 • p^ + nic) K"^ (7 " Pa + rUa) T'^] 

^D + E , (A26) 



where 



D = 



Qa 



Pa ■ q 



Pa,a {-g'^'q" + g'^q") 



X 



Tr 



dV 



(7 • Pc + rric) (7 • Pa + rua) T'^jt, + (7 • Pc + rrie) TsblpK 

OPa 



db 



+ (7 • Pc + rric) Tdb (7 • Pa + "^a) 



db 



dfa 



(A27) 



and 
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Qc 



Pc-q 



Tr 



^'dhlp^db (7 • Pa + ma) + (7 • Pc + m^) (7 • Pa + m„) T^^,, 

dpi 



+ (7 • Pc + "^c) Tdft (7 • Pa + rUa, 



db 



dp, 



Similarly, 



Tr [(7 • Pd + rud) Tea (7 • Pfe + "^fe) K'^c\ + Tr [(7 • p^ + rud) (7 • Pfe + ^b) T 



where 



and 



Now, 



where 



and 



Qb 



Pb ■ q 



-Pb. 



Tr 



(7 • Pd + rud) (7 • Pb + rUb) + (7 • Pd + rUd) T^alp^'cc 
dPb 

dV 



+ (7 • Pd + m^ ^ca {'y-Pb + rrib) 



dpt 



^ Qd ( ua u8 c 

G = Pd,a [9 q^ -9 ^q 



Pd-q 



X 



Tr 



^'calli^ca {l-Pb + mb) + (7 • Pd + rUd) (7 ■ Pb + rrib) T'^, 

or 



dp^d 



+ (7 • Pd + rUd) Tea (7 • Pb + rUb) 



dp^d 



J2 \Mof ^Ror S 



R = ^T: [(7 • Pd + md) Tea (7 • P6 + mfe) r'J , 



S = |Tr [(7 • pc + mc) Tdfe (7 • p„ + m„) r^J 
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Therefore, 



d\Mo\ 



rTr 



or 



(7 ■ + rUc) (7 ■ Pa + ma) T'^b + {l-Pc + rric) TdblaKb 
+ (7 ■ + rric) Tdb (7 ■ Pa + rua. 



or 



db 



(A35) 



d\Mo\ 
dPb 



dTca 

{l-Pd + rud) {l-Pb + mb) T'^a + {l-Pd + rrid) T^alaKc 

or 

+ {1 ■Pd + rrid) Tea {1 ■Pb + mb) ^-f 



(A36) 



dp"^ 



rTr 



dV 



T'dblaTdb (7 ■ Pa + ma) + (7 " Pc + m^) (7 " Pa + "^a) T'^b 

dV 

+ (7 ■ + mc) Tdb (7 ■ Pa + ma) 

OPc 



(A37) 



and 



d\Mo\ 
dPd 



rTr 



dTca 

T'calaTca {l ' Vd + md) + (7 " Pd + "^d) (7 " Pb + "^b) 



+ (7 ■ Prf + "^rf) Tea (7 ■ Pfe + "^6^ 



dT' 



(A38) 



Finally, 



'^il {Pi-(l){Pj-(l)^'''^ dp 



d\M,_ 



X 



(A39) 



where the sum in the second term runs over all fermions and Q'j = —Qj for the incoming 
fermions and Q'^ = Qj for the outgoing fermions. The first term will lead over the second 
one, in an expansion in powers of the virtual photon four-momentum. The above equation, 
together with 



E L,Ll = L 



(A40) 



proves Eq. @). 



25 



REFERENCES 

[1] See, for example, Proceedings of the Workshop on Dilepton Production in Relativis- 
tic Heavy Ion CoUisions, March 2-4, 1994, GSI, Darmstadt, Germany; edited by H. 
Bokemeyer. 

[2] G. Roche et al, Phys.Rcv.Lctt. 61, 1069 (1988);G. Roche et ai, Phys. Lett. B 226, 
228 (1989) ;C. Naudet et aL, Phys. Rev. Lett. 62, 2652 (1989). 

[3] L. Xiong, Z.G. Wu, CM. Ko and J.Q. Wu, Nucl. Phys. A 512, 772 (1990). 

[4] Gy. Wolf, G. Batko, W. Gassing, U. Mosel, K. Niita and M. Schafer, Nucl. Phys. A 
517, 615 (1990) . 

[5] Kevin Haghn and Gharles Gale, Phys. Rev. C 49, 401 (1994). 

[6] W.K. Wilson et al., Phys. Lett. B 316, 245 (1993). 

[7] Kevin Haghn, Joseph Kapusta, and Gharles Gale, Phys Lett. B224, 433 (1989). 

[8] M. Schafer, T. S. Biro, W. Gassing, and U. Mosel, Phys. Lett. B221, 1 (1989). 

[9] Kevin L. Haglin, Ann. Phys. (N. Y.) 212, 84 (1991). 

[10] L.A. Winckelmann, H. Stocker and W. Greiner, Phys. Lett. B 298, 22 (1993). 

[11] See, for example, J. D. Bjorken and S. D. Drell, Relativistic Quantum Mechanics 
(McGraw-Hill, New York, 1964). 

[12] R. Riickl, Phys. Lett. B 64, 39 (1976). 

[13] P. Lichard, Preprint SUNY-NTG-94-16 and Physical Review D, in press. 
[14] V. Balek, N. Pisiiova and J. Pisiit, Acta Phys. Slovaca 41, 224 (1991). 
[15] C. Gale and J. Kapusta, Phys. Rev. C 35,2107 (1987). 
[16] G. Gale and J. Kapusta, Phys Rev. C 40, 2397 (1989) . 

26 



[17] E. Byckling and K. Kajantie, Particle Kinematics (Wiley, New York, 1973). 

[18] L. Xiong, J. Q. Wu, Z. G. Wu, C. M. Ko, and J. H. Shi, Phys. Rev. C 41, R1355 (1990). 

[19] K. Haglin, C. Gale, and V. Emel'yanov, Phys. Rev. D 47, 973 (1993). 

[20] J. Cleymans, V. V. Goloviznin, and K. Redlich, Phys. Rev. D 47, 989 (1993). 

[21] R. Tabti, H. Eggers, and C. Gale, in preparation. 

[22] R.L. Hatch, Ph.D Thesis, Cahfornia Institute of Technology, 1980, (unpublished). 
[23] We have used DLS filter version 2.0, provided by the experimental collaboration. 
[24] A. I. Titov, B. Kampfer, and E. L. Bratkovskaya, JINR preprint E2-94-57, 1994. 

[25] E. L. Bratkovskaya, O. V. Teryaev, and V. D. Toneev, preprint DOE/ER/40561-167- 

1NT94-15-01, 1994. 



27 



FIGURES 



FIG. 1. The leading contributions for dilepton radiation in the reaction a + b^c + d + e^e . 

FIG. 2. The ratio R = — J^irtuai / ~ '^leai ^^'^ dependence on the center-of-mass dilepton 
momentum q = \q\ and hadron scattering angle 9 for (a) pn, and (b) pp processes at kinetic energy 
4.9 GeV and invariant mass M = 0.2 GeV. 

FIG. 3. The ratio of the differential cross section Ri = (^j^g^^) virtual / (^]g^^)reai as a function 
of the center-of-mass dilepton momentum q for (a) np, and (b) pp processes, at kinetic energy 4.9 
GeV and for dilepton invariant masses M= 0.01, 0.2, 0.6 and 0.8 GeV, respectively. 

FIG. 4. The differential cross section for dilepton production from (a) pn and (b) pp collisions 
versus the dilepton invariant mass M. The solid lines are results using the virtual photon formalism 
described in the text, the dotted lines are the results obtained using Riickl's formula together with 
the current for real photons. Results for kinetic energies -E'kin= 1-0, 3.0, and 4.9 GeV are shown. 
For clarity of presentation, the results have been multiplied by a scaling factor. The scaling factor 
is 1 at 1 GeV, 10 at 3.0 GeV and 100 at 4.9 GeV. 

FIG. 5. The same as in Fig. 4, but with the correction of the acceptance filter. The solid 
symbols are with the formalism of Eq. (^) together with the current of Eq. (^), while the open 
symbols are the results of using Riickl's formula together with the current for real photons. 

FIG. 6. The ratio of dilepton production cross sections in pp and np reactions, 
^ ~ (^)pp/(^)pn' ^ function of invariant mass M at kinetic energies 1.0, 2.0, 3.0 and 
4.9 GeV, from bottom to top. 
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FIG. 7. We plot the differential cross section for production of lepton pairs of invariant M in 
the cases oi pn (left panel) and pp collisions (right panel) at incident kinetic energies of 1, 3, and 
4.9 GeV (bottom to top curves, respectively). The full curves represent calculations done with Eq. 
([T^ , the short-dashed curves are generated with the formula of Riickl, together with the current 
for real photons, and the long-dashed curves are done using Eq. (^) together with the current of 
Eq. (|^. The scaling factors are as in Fig. |^. 

FIG. 8. The same as Fig. 0, but with the DLS acceptance corrections. The open triangles are 
the results of using Riickl's formula together with the current for real photons. The open circles 
are associated with Eq. (|8|) together with the current of Eq. (^ . The solid triangles are generated 
withEq. M). 



FIG. 9. Same caption as Fig. ^, only this time the complete phase space approach of section 



IV is used. 



FIG. 10. We plot the polar anisotropy coefficient, as defined in the text, against lepton pair 
invariant mass for incident kinetic energies of 1, 2, 3, and 5 GeV in pn collisions. We use the 
current of Eq. (^ . 



FIG. 11. Same caption as Fig. 10. We use here the current of Eq 



FIG. 12. Same caption as Fig. but for pp collisions. 
FIG. 13. Same caption as Fig. but for pp collisions. 
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